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ABSTRACT 

We investigate the properties of 1088 Lyman Break Galaxies (LBGs) at z ~ 3 selected from a ~ 2.63 deg 2 
sub-region of the First Look Survey field using the ground-based multi-color data and the Spitzer Space Tele- 
scope mid-infrared data at 3-8 and 24 /im. With the wide area and the broad wavelength coverage, we sample a 
large number of "rare" M-band dropouts which are massive (M» > 10 1 1 Mq), allowing us to perform a statistical 
analysis of these subsets of LBGs that have not been studied in detail. Optically bright (Rab < 24.5 mag) LBGs 
detected in mid-infrared {Si^^m > 6/xJy) reside at the most massive and dusty end of the LBG population, with 
relatively high and tight M/L in rest-frame near-infrared. Most infrared-luminous LBGs {Su^m > 100/iJy) are 
dusty star-forming galaxies with star formation rates of 1 00-1 000 M Q /yr, total infrared luminosity of > 10 12 Lq. 
By constructing the UV luminosity function of massive LBGs, we estimate that the lower limit for the star for- 
mation rate density from LBGs more massive than 10 n M Q at z ~ 3 is > 3.3 x IQ^MQ/yr/Mpc 1 , showing 
for the first time that the UV-bright population of massive galaxies alone contributes significantly to the global 
star formation rate density at z ~ 3. When combined with the star formation rate densities at z < 2, our result 
reveals a steady increase in the contribution of massive galaxies to the global star formation from z = to z ~ 3, 
providing strong support to the downsizing of galaxy formation. 

Subject headings: cosmology: observations - galaxies: evolution - galaxies: high-redshift - galaxies: stellar 
content - galaxies: starburst - infrared: galaxies 



1. INTRODUCTION 

Stellar population analysis of local galaxies shows that stars 
in massive galaxies formed early within a short period in the 
history of the universe, while low-mass galaxies went through 
late, slow star formation (Thomas et al. 2005; Heavens et 
al. 2004; Panter et al. 2006). The appearance of passively 
evolving early-type galaxies above z > 1 (Cimatti et al. 2002; 
Forster-Schreiber et al. 2004) supports the idea of early star 
formation in massive systems. The evolution of specific star 
formation rate, i.e., the star formation rate per unit stellar 
mass shows that the star formation in most massive galax- 
ies (M* > 10 11 M a ) have nearly completed by z ~ 1.5 while 
that for less massive galaxies continue to date (Papovich et al. 
2006). These observations support the "downsizing" scenario 
of galaxy formation, commonly expressed as the decrease of 
stellar masses of galaxies in which vigorous star formation 
occurs with decreasing redshift (Cowie et al. 1996; Kodama 
et al. 2004; Juneau et al. 2005). 

Attempts have been made to explain this seemingly anti- 
hierarchical behavior of galaxy formation using hierarchical 
galaxy formation models. In an attempt to explain downsiz- 
ing within the frame of semi-analytic model, Neistein, van 
den Bosch, and Dekel (2006) distinguishes "archaeological 
downsizing "(ADS) versus "downsizing in time"(DST). In 
ADS, the downsizing is characterized as the built-up of stellar 
masses where massive galaxies are assembled through merg- 
ers of less massive galaxies which formed stars early. In com- 
parison, DST is characterized by the formation and assem- 
bly of massive systems at early epoch. Semi-analytical mod- 
els seem to be able to explain ADS naturally, while compli- 
cated baryonic processes are necessary to explain DST. Re- 
cent semi-analytic models have shown that the inclusion of 
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AGN feedback effect successfully reproduce the observed de- 
crease of actively star-forming system, yet more understand- 
ing about the observed quantities and cooling processes are 
needed. On the other hand, some other models do not seem 
to have much difficulty explaining the abundance of massive 
galaxies at high redshift (e.g., Nagamine et al. 2005). 

Recent observational evidences point toward the DST-type 
of downsizing which suggests that there should be "massive", 
star-forming galaxies at high redshift. Some of the red, mas- 
sive galaxies are found to be dusty star-forming galaxies, and 
submm galaxies are thought to be massive galaxies in for- 
mation at high redshift (Smail et al. 2002). Still, it is not 
clear if they account for the whole star formation activity of 
massive galaxies. To complete the picture of the downsizing 
galaxy formation and to provide additional constraints to re- 
fine the models, we need to identify and study other actively 
star-forming massive galaxies, such as Lyman break galaxies 
(LBGs). 

Lyman break galaxies, galaxies selected by the continuum 
break at the Lyman limit, are the most commonly studied 
high-redshift star-forming galaxies. Since their selection tech- 
nique requires only optical imaging observation, many LBGs 
up to z ~ 6 have been discovered and studied to date (Stei- 
del et al. 1999; Papovich et al. 2001; Shapley et al. 2001; 
Giavalisco 2002; Ouchi et al. 2004; Bouwens et al. 2004). 
These studies reveal that typical z ~ 3 LBGs are galaxies with 
large ongoing star formation of 10-100 M Q /yr (Papovich et 
al. 2001; Shapley et al. 2001). Typical stellar masses of 
z ~ 3 LBGs are found to be mostly of order of 10 10 M Q or 
less from the spectral energy distribution fitting using optical 
and near-infrared photometry, although LBGs with masses up 
to 10 11 M Q were also found (Papovich et al. 2001; Shap- 
ley et al. 2001). Some simulations suggest that nearly 50% 
of > lO n M galaxies could be detected using the rest-frame 
UV selection criteria (Nagamine et al. 2005), but the observed 
fraction of LBGs among massive galaxies is less than 20% 
(van Dokkum et al. 2006). 
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Recently, the addition of Spitzer mid-infrared (MIR) ob- 
servations not only reduced the uncertainties in derivation of 
stellar mass (Shapley et al. 2005), but also enabled the estima- 
tion of the infrared luminosities of LBGs. Huang et al.(2005) 
have defined Infrared Luminous LBGs (ILLBGs) as LBGs de- 
tected in the Spitzer MIPS 24 pm with fu^m > 60^Jy. The 
infrared luminosities of ILLBGs are estimated to be larger 
than 10 12 L Q , and they contribute ~ 5% of total LBG pop- 
ulation. The stellar masses of ILLBGs exceed 5 x 10 10 M Q 
(Rigopoulou et al. 2006), and their star formation rates in- 
ferred from infrared luminosity are as high as ~ 1000M Q /yr 
that is sufficient to evolve into present-day giant ellipticals. 
Several massive LBGs, not necessarily ILLBGs, are also iden- 
tified in a part of the Extended Groth Strip covering ~ 227 
arcmin 2 (Rigopoulou et al. 2006). These massive LBGs are 
found to have a star formation rate (SFR) of > 100M Q /yr. 
These results suggest that ILLBGs and massive LBGs are im- 
portant indicators of the star formation activity in massive 
galaxies at high redshift, along with the distant red galaxies 
(DRGs; Franx et al. 2003) and the submm galaxies. 

However, the number of z ~ 3 massive LBGs and ILLBGs 
is still small due to the difficulty of covering large areas with 
expensive M-band observations in order to discover such rare 
objects. For example, only 5 massive (> 10 n Mo) LBGs are 
found in Rigopoulou et al.(2006). At z ~ 3, the largest LBG 
sample (2347 photometric LBGs) until now is gathered from 
1 1 separate fields of ~ 0.38 deg 2 in total (Steidel et al. 2003). 
LBGs selected in the ESO Deep Public Survey (Hildebrandt et 
al. 2007) covers as much as ~ 1.75 deg 2 at bright end (30' x 
30' x 7; R AB < 24.0), but until now, the survey consists of 
optical data only. 

In this study, we enlarge the z ~ 3 massive LBG/ILLBG 
sample and investigate their properties using Spitzer First 
Look Survey. The survey is not as deep as other multi- 
wavelength surveys, however, it is wide enough (~ 4.3 deg 2 , 
of which 2.63 deg 2 was used in this study) to collect the rare 
bright LBGs. LBGs selected in optical bands (with /?-band 
magnitude cut of Rab < 24.5 mag) are combined with detec- 
tions in Spitzer IRAC (3.6-8/xm) and MIPS (24/xm). By com- 
paring observed flux with model galaxy spectra, we estimate 
ages, the amount of dust extinction, and stellar masses of the 
LBGs. We also present 2\[im detected LBGs in our study, 
and compare them with the ILLBGs in Huang et al.(2005). 

The data set we used is described in §2. We describe our 
selection criteria for z ~ 3 LBG candidates in §3. The anal- 
ysis on the stellar population of LBGs follows in §4, and the 
key results are presented in §5, such as the luminosity func- 
tion, the stellar mass properties, and dust properties. We also 
discuss the implication of our study on the galaxy evolution, 
in particular regarding the "downsizing" picture of galaxy for- 
mation in §5. Finally in §6, the conclusion follows. Through- 
out this paper, we use a cosmology with VLm = 0.3, ft\ = 0.7, 
and Hq = lOkm/s/Mpc. All the magnitudes are specified in 
AB system, unless noted. 

2. DATA 

The LBG candidates are selected from the extragalactic 
component of the Spitzer First Look Survey (XFLS). XFLS 
comprises of InfraRed Array Camera (IRAC, 3.6, 4.5, 5.8, 
8.0^m; Fazio et al. 2004) and Multiband Imaging Photome- 
ter for Spitzer (MIPS, 24, 70, 160^; Rieke et al. 2004) 
imaging observations over the ~ 4.3 deg 2 field centered at 
RA= 17H8"W, DEC= +59°30'00". With an effective expo- 
sure time of 1 minute per pixel, the main survey field (~ 4.3 
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FIG. 1 . — The area coverage maps of various datasets in XFLS. The small- 
est region marked with thickdashed line in the center represents verification 
field for IRAC (~ 900 arcmin 2 ) that is slightly deeper than other parts of 
XFLS. The skewed rectangle (dash-dot-dash line) is the MIPS verification 
strip. The "central 1 deg 2 " region is drawn as a thick square in the center. 
The shaded area outside the central 1 deg 2 region represents the outer field 
(~ 1.63 deg 2 ) where bright w-dropouts are selected. The IRAC main field 
of XFLS is specified with the thin dashed line at outermost. Thin Solid line 
shows our «-band coverage. 

deg 2 ) has 5 a flux limits of 20, 25, 100, and lOO^Jy at wave- 
lengths of IRAC 3.6, 4.5, 5.8, and 8.0/xm respectively. The 
flux limit of MIPS 24/xm is about ~ 300/Jy (5 a). The cen- 
tral part of the FLS (~ 900 arcmin 2 ) was observed deeper, 
with 10-minute integration times per pixel. This "verifica- 
tion" field has sensitivities of 10, 10, 30, and 30/iJy at IRAC 
wavelengths (5 a; Lacy et al. 2005). For the MIPS verifi- 
cation field, 3 a flux limits are 90y(iJy at 24^m (Fadda et al. 
2006), 9 and 60 mJy at 70 and 160^m respectively (Frayer et 
al. 2006). 

In order to select M-dropout objects in the XFLS, we used 
u, g, and /?-band ground-based images acquired by three dif- 
ferent wide-field instruments. We obtained deep M*-band im- 
ages with MegaCam on the Canada-France-Hawaii Telescope 
(CFHT) 3.6-m telescope (Shim et al. 2006). CFHT M *-band 
filters are slightly redder than popular SDSS M-band. There- 
fore using the similar color selection criteria, we expect that 
the M*-band dropouts are biased towards higher-redshift ob- 
jects than previously selected LBGs. The central 1 deg 2 field 
was observed to the depth of u* ^26.2 mag (AB, 5(7 in 3" 
diameter aperture). In addition to the central part, the whole 
FLS area was covered by M*-band imaging with a shallower 
depth (~24.5 mag). Detailed information about the observa- 
tions, photometry and the properties of the dataset are pre- 
sented in a separate paper (Shim et al. 2006). Deep g-band 
image was obtained together with the «*-band observation 
over the central 1 deg 2 (g <~ 26.5 mag). 

The g'-band images were obtained over ~ 2.63 deg 2 , in- 
cluding the central 1 deg 2 area, using Large Format Cam- 
era (LFC; Simcoe et al. 2000) at the Palomar 5-m tele- 
scope (Storrie-Lombardi et al., in preparation). The LFC 
g'-band images reach to the depth of 24.5 mag. For the re- 
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gion surrounding the central 1 deg 2 area, we use the LFC 
data for g-band photometry. The u*,g catalogs were ex- 
tracted through dual-mode photometry of SExtractor (Bertin 
& Arnouts 1996), using the g-band image as a reference im- 
age. The 7?-band catalog of the FLS (Fadda et al. 2004) were 
matched with u*, g-band catalogs using the matching radius 
of 0.8", to obtain the tf-band flux. The depth of the tf-band 
image varies depending on the field location. Still, it reaches 
R(Vega) ~ 24.5 mag at the deepest central field. Our method 
to select z ~ 3 LBGs with optical colors is described in detail 
in Section 3. 

Other ancillary datasets used to investigate the properties 
of M-dropout objects include i', J, and /T,-band images. LFC 
/'-band images cover the same area with the LFC g'-band im- 
ages with the depth of ~ 24 mag. For the central part less than 
1 deg 2 , we also used NIR (J(Vega) ~ 21 mag, K s (Vega) ~ 20- 
21 mag) data, obtained with WIRC on the Palomar 5-m (K s - 
band), and FLAMINGOS on the KPNO 4-m (/-band; Choi 
et al., in preparation). The whole FLS field was covered by 
1.4 GHz radio observations with VLA (flux limit at 90 fiiy, 
Condon et al. 2003). In this study, spectroscopic redshifts ac- 
quired with DEIMOS/Keck for -1300 objects at z - 1 (Choi 
et al. 2006) were used to check if there is any low-redshift 
interlopers within the LBG sample. 

3. LBG SAMPLE SELECTION 
3.1. Photometric Selection Criteria for u-dropouts 

We select M-dropout objects from the 2.63 deg 2 area that is 
covered by u*, g, R, and /-bands. This area is illustrated as 
a shaded region in Figure 1. We define M-dropouts as objects 
that satisfy /?-band magnitude cut and color-cut criteria in g — 
R versus u* -g diagram. The method is in principle identical 
to the popular "drop-out" method for selecting LBGs (e.g., 
Steidel & Hamilton, 1993). 

The adopted selection criteria can be summarized as fol- 
lows: 

u-g> 1.4, 
u-g>3(g-R)-Q.U, 
g-R< 1.2, 

R < 24.5 (for central 1 deg 2 ), R < 23.5 (for the outer area) 

Different /?-band magnitude cuts are adopted taking into ac- 
count the difference in image depth (mainly g'-band) between 
the central and the outer regions. For the central 1 deg 2 field, 
we apply a magnitude cut of R < 24.5 magnitude. For the 
outer 1.63 deg 2 area surrounding the central region, we adopt 
R < 23.5 magnitude cut. Note that u, g, and R magnitudes are 
corrected for the galactic extinction (Shim et al. 2006; Fadda 
et al. 2004). 

The colors of star-forming galaxies at 2.8 < z < 3.6 falls 
into the selection box (Figure 2). The lines in the upper left of 
the color-color plot represent the location of model galaxies 
with constant star formation rate at z ~ 3, of various ages and 
extinction values (see Figure 2 caption). Note that the lines 
represent galaxies with the solar metallicity, and the metallic- 
ity lower than the solar would move the lines to the left. 

The color selection criteria can select star-forming galax- 
ies with constant star formation rate at the age of — 10 Myr 
to — 2 Gyr. The criteria can also select star-forming galaxies 
with exponentially decaying star formation rate at the age up 
to the exponential time scale r. Regardless of its star forma- 
tion history, a galaxy younger than — 100 Myr falls into the 
selection box. When there is internal extinction, the limit on 
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Fig. 2. — The g-R versus u* —g color-color diagram of galaxies in 
XFLS. Using the expected tracks of galaxies at the redshift range of 
2.8 < z < 3.6 with different age/star formation history/dust extinc- 
tion, we identify the boxed region in the upper left as the place where 
the LBGs reside. The solid lines indicate galaxies with constant star 
formation, and the dashed lines indicate galaxies with single-burst 
star formation history. Different line colors indicate different ages, 
blue for 50 Myr-old galaxies and red for 1 Gyr-old galaxies. The 
shifts in x-axis direction represent the amount of dust extinction. 
Galaxies with constant star formation younger than 2 Gyr can be 
selected using the specified criteria. Young galaxies that have burst- 
like star formation history can also be selected, but old galaxies at 
their passively evolving stages (red, dashed line on the right) cannot 
be selected with this criteria. 

age of galaxies that can be selected with the criteria changes. 
For example, the star-forming galaxies with age less than 500 
Myr can be selected if E(B-V) < 0.3. Passively evolving 
galaxies as old as — 1 Gyr at z ~ 3, which are growing old 
after a single burst, cannot be selected by this criteria (red, 
dashed line on the right in Figure 2). 

After selecting M-dropouts from the color-color space, we 
remove spurious objects through the visual inspection. We 
find that spurious objects are mostly lying near the edge of 
the image, or too close to bright stars. On the other hand, 
we exclude objects that are likely to be a low-redshift inter- 
lopers based on either spectroscopic redshifts or photometric 
redshifts (see SED fitting section on Section 4.1). 

Finally as a result, we identify 1088 M-dropouts in our 
search area (925 objects with R < 24.5 mag in the central 1 
deg 2 , and 163 objects with R < 23.5 mag in the outer 1.63 
deg 2 ). Table 1 shows the result. As mentioned in the above 
paragraph, the numbers in Table 1 are corrected for possible 
low-redshift interlopers based on their photometric redshifts 
or spectroscopic redshifts. At R > 23 mag, the number density 
of M-dropouts of our study is consistent with the LBG surface 
density from previous studies (e.g., Steidel et al. 1999). At 
22.5 < R < 23, our surface density appears to be about twice 
higher than those found in Steidel et al.(1999), although 1-<t 
error bars of the two numbers overlap. On the other hand, our 
number is about three times lower than those found by Hilde- 
brandt et al. (2005). The possible explanations for the dis- 
crepancy are (i) low-redshift interlopers; (ii) the cosmic vari- 
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TABLE 1 

Surface Density of w-band Dropout Objects 



Rab magnitude 


N " 


number density b 


Steidel+99 c 


Hildebrandt+05 d 






(deg- 2 ) 


(deg- 2 ) 


(deg- 2 ) 




22.0-22.5 


11 


4.18±1.26 (1.72) 






22.5-23.0 


35 


13.31±2.25 (4.33) 


6.77±3.38 


34.2±10.8 


23.0-23.5 


188 


71.48±5.21 (20.6) 


71.06±13.54 


72±18 


23.5-24.0 


247 


247±15.7 (35.5) 


294.4±37.22 


234±28.8 


24.0-24.5 


607 


607±24.6 (81.9) 


656.50±54.14 


540±54 


1088 



a w-band dropouts are selected in the effective area of 2.63 deg 2 when R < 23.5 mag, and 1 deg 2 when R < 24.5 mag. 

b The errors are given for the Poisson noise only, and in the parenthesis, the errors considering both the Poisson noise and the galaxy clustering 
noise (a 2 = oj> oisson + o 1 chaering ) are given. To calculate the clustering error, we used the approximation formula in Peebles (1975), adopting the 
amplitude of the angular two-point correlation function of Iab < 24.5 LBGs in Foucaud et al.(2003). The value in Foucaud et al.(2003) was 
consistent with our own derivation within ~ 30 %. 



The surface number density of z ~ 3 LBGs from Steidel et al.(1999). The errors given here represents the Poisson error. Steidel et al.(1999) 
used TZ filter, which is slightly different from R filter with TZab -Rab = -0.11 (Foucaud et al. 2003). 

d The surface number density of z ~ 3 LBGs (w-dropouts) in CDF-South (Hildebrandt et al. 2005). The areas used in their studies are ~ 900 
arcmin 2 . 



TABLE 2 
Subsets of LBG Sample 



Group total 3.6/im 4.5/im 5.8/im 8.0/im 24/im 



ver (1) (R <24.5) 189 22 21 6 7 4 
cent(2) (R <24.5) 736 21 20 2 3 2 
outer (3) (R <23.5) 163 20 18 4 4 6_ 

NOTE. — (1) The IRAC verification field (~ 0.22 deg 2 ) (2) The central part excluding the verification (i.e., effective area: ~ 0.77 deg 2 ) (3) 
The outer ~ 1.63 deg 2 field 



ance; and (iii) AGN contribution (see Section 5.1 and Figure 
6). Spectroscopic observation of bright LBGs should be able 
to determine which explanation is right. Our study has an 
advantage over previous ones in terms of the area coverage 
(~ 2.63 deg 2 ), thus the Poisson errors in the number density 
estimates are significantly reduced. 

3.2. Mid-Infrared Detection of LBGs 

The infrared fluxes of M-dropout objects are obtained from 
the published IRAC/MIPS source catalogs when available 
(Lacy et al. 2005; Fadda et al. 2006). If the object flux is 
not available, we measured the object flux through empirical 
PSF fitting, applying aperture correction values from Lacy 
et al.(2005). To reduce uncertainties in flux measurement, 
we use only objects whose fluxes exceed 3-a flux limits in 
IRAC. The number of M-dropouts detected in each IRAC band 
is summarized in Table 2. Here, we define LBGs detected 
over 3-a limit in 3.6 /mi image as an IRAC-detected LBGs 
(hereafter, IRAC LBGs) since all the objects detected in other 
IRAC bands are also detected in 3.6 /nm. 

The number of IRAC LBGs is dependent on the depth of the 
Spitzer image. For example, the number of detections in Table 
2 is relatively small compared to the deeper surveys because 
of the shallower depth of our survey (e.g., Barmby et al. 2004; 
Huang et al. 2005). Note that, at the same infrared flux limit, 
the number of IRAC detections in our study is consistent with 
the other studies. For example, we find 22 IRAC LBGs above 
the 3 a flux limit of 6 /iJy (3.6/im) in the FLS verification 
strip (900 arcmin 2 ). On the other hand, Huang et al.(2005) 
have identified 8 LBGs and 3 QSO/AGNs in 3.6 /im image at 
the same flux limit, over 227 arcmin 2 (see Figure 1 in Huang 
et al. 2005). If we scale the area coverage, we would expect 
that we could find ~ 6 objects in 227 arcmin 2 of Huang et 



al.(2005). These numbers are relatively consistent with each 
other. Huang et al.(2005)'s sample includes LBGs as faint 
as 24.5 < R < 25.5 while we do not. This /?-band difference 
may cause a slight discrepancy between the two numbers, but 
the difference is almost negligible for the LBGs bright in the 
rest-frame NIR wavelengths which also turn out to be mas- 
sive LBGs. Given the bright rest-frame NIR flux, we find that 
IRAC LBGs belong to massive, IR-bright end of LBGs (see 
Section 5.2 for more detail). Our sample is biased against 
UV-faint LBGs due to the bright /?-band magnitude cut, but 
the above comparison suggests that the bias does not affect 
our analysis of massive LBGs since they are mostly bright in 
IRAC wavelengths. 

For the MIPS detection, we matched our IRAC LBGs with 
MIPS 24 fim sources in Fadda et al.(2006). In the verification 
field, 6 IRAC LBGs found matches in the 24 /im catalog (over 
5 a detection, corresponding to S24 ~ 150/iJy). From the shal- 
lower parts of the FLS field, 6 objects are found to be in the 
24 fim catalog (5 a, 300/dy). We call these LBGs detected 
in 24 fim as "24 fim LBGs" 4 . In total, we have twelve 24 
^m LBGs from 1088 M-dropouts in our sample. Similarly to 
the case of IRAC LBGs, MIPS -detection is dependent on the 
depth of the image, and the number of 24 /im detections can 
certainly go up if we had a deeper MIPS image. The 24 /im 
LBGs are IR-luminous sub-population of z ~ 3 IRAC LBGs 
with L IR > 10 12 L Q (See section 5.3 for more detail). 

4. SED FITTING 

We have performed SED fitting in order to derive vari- 
ous parameters, especially focusing on the stellar masses of 

4 Among non-IRAC LBGs, 2 objects found their matches in 24 fim im- 
ages. However, we didn't include them in 24 fim LBG sample since their 
properties are hard to constrain with limited optical photometry. 
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FIG. 3. — The SED fitting results for a representative subset of our LBG sample. The x-axis is the observed wavelength, and the photometric data points are 
indicated for u, g, R, i, J, K s , IRAC channel 1, 2, 3, and 4 at 0.374, 0.487, 0.651, 0.768, 1.26, 2.16, 3.6, 4.5, 5.8, and 8.0 fim. The object IDs are drawn from 
i?-band catalogs of Fadda et al.(2004). The solid line is the best-fit template using all the photometry data-points. For some objects, a dashed line that indicates 
the best-fit template without MIR data-points is also overplotted. 



LBGs. The SED fitting method employs the same method- 
ology as previous studies (e.g., Papovich et al. 2001; Shap- 
ley et al. 2001). Various galaxy SED templates were gen- 
erated using a stellar population synthesis model (Bruzual & 
Chariot 2003), and fitted to the observed SED. We used all 
the available photometric data points from w-band to IRAC 8 
/zm. Note that if we have u,g, and /?-band data only, we did 
not perform the SED fitting since it was very difficult to gain 
meaningful fit-results from three photometric data points. For 
optical fluxes, we used MAG_AUTO in SExtractor, while in- 
frared fluxes were measured from the PSF fitting or drawn 
from the catalog (Lacy et al. 2005) as described in the pre- 
vious section. In the NIR, the LBG candidates were matched 
with J- and tfj-band catalogs (MAG_AUTO from SExtractor; 



Choi et al. 2007, in preparation) within the matching radius 
of 1.2". 

For the generation of various galaxy SEDs, we limited our 
templates to a specific star formation history. The galaxies 
have the exponentially decaying form of (0(f) oc expi-t/r)) 
star formation with r=10Myr, lOOMyr, 300Myr or a constant 
star formation rate. In our case, the best-fit SED of almost 
all the LBGs gives a constant star formation template (cf. 
Rigopoulou et al. 2006). The metallicity and the initial mass 
function are fixed to certain values/form. The adopted metal- 
licities are either 0.2 Z Q or 1 Z©. The IMF was fixed as the 
Salpeter IMF between 0.1 M Q and 100 M Q . The stellar pop- 
ulation was evolved using Padova 1994 stellar evolutionary 
tracks. The ages of a galaxy, f, was chosen from ^22 possi- 
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FIG. 4. — Le/h The photometric redshift distribution of our LBG sample derived from SED fitting. Dotted line is an approximation of the distribution as a 
Gaussian function. Mean redshift (z) is ~ 3.2, and the standard deviation of the distribution is ~ 0.14. Right: The comparison between spectroscopic redshifts 
and photometric redshifts of LBGs in the Westphal field (Steidel et al. 2003). Open circles represent 8 micron detected LBGs from Rigopoulou et al. (2006), 
while filled circles indicate other LBGs with at least one NIR (J,K S ) photometry. Among the 8 micron LBG sample, only 6 objects have NIR photometry (points 
with both filled and open circles). 



ble values between 10 Myr to the age of the universe at the 
corresponding redshift. The galactic reddening of the gener- 
ated spectra is following the dust extinction law for starburst 
galaxy (Calzetti et al. 2000). We also added the flux suppres- 
sion shortward of Ly a forest due to the intergalactic medium 
following Madau et al.(1995). 

The best-fit parameters - photometric redshift z, extinction 
parameter E(B - V), stellar mass, age, star formation his- 
tory and metallicity - are obtained by minimizing the error- 
weighted x 2 value expressed as x 2 = X); filter ^" b ~i F "^ . Note 
that ages, star formation histories, and metallicities are chosen 
from discrete values as explained above. 

An illustration of SED fitting for some galaxies is presented 
in Figure 3, including SEDs of possible low-redshift interlop- 
ers. 

4.1. Photometric Redshift 

The distribution of photometric redshifts of M-dropouts is 
presented in Figure 4a. The median redshift of LBGs is 
(z) ~ 3.2, and the standard deviation of the photometric red- 
shift distribution of LBGs is ~ 0.14. The M-dropouts are dis- 
tributed mostly at 2.5 < z < 3.8, but there is a low redshift tail 
to the distribution. We consider those having low photometric 
redshifts to be interlopers. 

For example, 5 objects out of the "initial" IRAC LBGs are 
turned out to have probable redshifts much less than 3. The 
best-fit SED results for these objects are heavily attenuated, 
moderately old galaxies at z ~ 1 (see the first object in the last 
row of Figure 3 as an example). Among these, 3 objects were 
identified to lie at z = 0.46,0.96, and 1.26 in match with the 
spectroscopic sample of z ~ 1 galaxies (Choi et al. 2006). The 
"possible contaminants" are characterized by red g-R and 
R-i colors, bright NIR flux (/), and decreasing MIR flux from 
3.6 fim to 8.0 which suggests that 1.6 /im H~ bump has not 
redshifted out of 3.6 fim. For LBG candidates with only g,R, 
and ;' detections, low-redshift interlopers are also identified as 
objects with red g-R and R- i colors. When there are even no 
i detection, it was difficult to determine which M-dropouts are 
low-redshift interlopers. However, the non-detection in /-band 



suggests that they have blue g-R colors which is indicative 
of LBGs atz~3. 

Through these breakdown of the low-redshift interlopers, 
we removed M-dropouts that have most probable photometric 
redshifts of z < 2. Our investigation shows that the fraction 
of possible low-redshift interlopers is higher at the brighter R- 
band magnitude bin. Among M-dropouts, the fraction of possi- 
ble interlopers with low photometric redshifts is ~ 25% at the 
brightest (22 < R < 22.5) bin, ~ 13% at the second brightest 
(22.5 < R < 23.0) bin. For the other bins, the fraction is of or- 
der of ~ 5%. Three spectroscopically confirmed low-redshift 
galaxies mentioned above and galaxies with low photometric 
redshifts are excluded in Table 1 and in the statistical analysis 
hereafter. A future spectroscopic study of bright LBGs will 
be able to determine how reliable the estimates of the number 
of interlopers are. 

In order to check the reliability of our photometric redshift, 
we tested our method on spectroscopically confirmed LBGs 
in the extended Westphal-Groth strip (Steidel et al. 2003; 
hereafter Westphal LBGs). For the photometric data, we used 
the optical, NIR, and MIR photometry of Westphal LBGs in 
Shapley et al.(2001) and Rigopoulou et al.(2006). All the 
LBGs used in this test have at least one photometric data point 
at wavelength longward of the optical (i.e., NIR or MIR). 

Figure 4b shows the comparison of our photometric red- 
shifts versus spectroscopic redshifts. Although there exist a 
few outliers, in overall the plot shows that the derived photo- 
metric redshift is quite reliable to about Az/z ~ 0.1, regard- 
less of whether we have NIR or MIR data. 

4.2. Stellar Mass 

Stellar mass is known to be a robust parameter that can be 
constrained relatively easily compared to other parameters, by 
being insensitive to the assumed star formation history (e.g., 
Papovich et al. 2001; Rigopoulou et al. 2006). Other con- 
straints, such as metallicity, are found to affect the derived 
stellar mass within a factor of 2-5 (Papovich et al. 2001). In 
addition to that, it has been reported that the inclusion of MIR 
photometry data points reduces the stellar mass uncertainties 
by a factor of 1.5-2 (Shapley et al. 2005). Here, we inves- 
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Therefore, we converted /?-band magnitude to the absolute 
magnitude with the K-correction consisting of the bandwidth 
dilation term only, and used it as the UV absolute magnitude. 
The UV luminosity derived this way samples different rest- 
frame UV wavelengths at different redshift, but the effect due 
to this should be less than a few tenths of magnitude. 

To derive the luminosity function, we used the y— method 
(Schmidt 1968; Lilly et al. 1995; Im et al. 2002). The number 
density of LBGs in each magnitude bin is calculated as 

<j)(M) [mag' 1 Mpc~ 3 ] = x ^ — !— 

while AM is the size of the magnitude bin. The maximum 
comoving volume, V mas is calculated by the equation 

/■mmfe,***) d y 

V nmx = I —dz 



max(zi,z m m) 



dz 



FIG. 5. — Left: The comparison of the stellar masses derived using spectroscopic redshifts and our method against that of Rigopoulou et al.(2()06). Right: The 
comparison of the stellar masses using photometric redshifts and that of Rigopoulou et al.(2006). The stellar masses estimated by two different methods agree 
with each other within the errors. 

tigate how our stellar masses fare with the values from other 
methods, and also how much the lack of spectroscopic red- 
shifts influences the stellar mass derivation. Again, we use 
the Westphal LBGs for this purpose (Rigopoulou et al. 2006). 

To inspect the uncertainty arising from photometric red- 
shift, we performed this comparison using spectroscopic red- 
shifts first (Figure 5 a). Then, we estimated stellar masses 
using photometric redshifts instead of spectroscopic redshifts 
(Figure 5b), and studied how the use of the photometric red- 
shifts affects the result. 

Figure 5a shows that our stellar masses derived with spec- 
troscopic redshifts are consistent with those from Rigopoulou 
et al. (2006) within < 0.2 dex. Therefore, we consider the 
inherent amount of stellar mass uncertainty due to the fit- 
ting method to be ~ 0.2 dex. Figure 5b shows that stellar 
masses derived with photometric redshifts are within 0.1-0.2 
dex from those derived with spectroscopic redshifts. Consid- 
ering that the uncertainty of stellar masses is dependent on 
the fitting method by about ~ 0.2 dex, the stellar masses de- 
rived with photometric redshifts and spectroscopic redshifts 
are consistent within the error from the fitting method. Two 
objects show relatively large discrepant values (objects at 
around [10.0, 10.5]), but this can be explained with a fail- 
ure of the fitting since these objects have large x 2 values. In 
addition, we have tested two-component fitting for several ob- 
jects (see last 2 objects in Figure 3) to see how the presence of 
an underlying old stellar population affects the derived stellar 
mass. For the two components, we used an old component 
(500Myr or lGyr passively evolving population after single 
burst) and a young component (constant star forming popu- 
lation younger than lOOMyr). The result shows <~ 20% of 
difference in the derived stellar masses compared to those de- 
rived using a single-component fitting. Since the discrepancy 
is relatively small, we used single-component fitting only. 

5. RESULTS 

5.1. Bright End of the Luminosity Function 

We have constructed the UV luminosity function of our 
LBGs using the entire sample corrected for interlopers (Ta- 
ble 1). At 2.5 < z < 3.8 where our LBGs are distributed, the 
central wavelength of /?-band corresponds to 1360A-1860A. 



where zi and zi indicate the lower and the upper limit of the 
redshift distribution of the galaxies. The z m in,z m ax values rep- 
resent the minimum and maximum redshift where the galaxy 
can be detected. dV /dz is the differential comoving volume 
at redshift z. 

Our result is consistent with previous results (Steidel et al. 
1999) at the magnitude range of -22.5 < M uv < -21.0 (Fig- 
ure 6). At the bright end, there are two noticeable results: (i) 
the decrease of error-bars, (ii) a clear excess of bright LBGs 
compared to the expected number from the best-fit Schechter 
function. First is due to the large area coverage of this study, 
since the error-bars in Figure 6 represent the Poisson error 
only. When the galaxy clustering is taken into account, the er- 
ror bar increases by up to a factor of 2 over the Poisson statis- 
tics. Even after the clustering effects are taken into account 
(Peebles 1975; please refer to the clustering error in Table 1), 
our error bars are smaller than those of previous studies. 

The excess of bright LBGs are quite interesting, although 
it has been expected from the high surface density of LBGs 
at the bright /?-band magnitude bin (Table 1). The apparent 
/?-band magnitudes of LBGs with M uv < -23.0 mag are all 
between 22 < R < 23. Until now, few LBGs have been dis- 
covered over M uv = -23.0 mag, including six "very" bright 
(M uv <~ -25.4 mag) M-dropouts discovered in SDSS DR1 
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FIG. 6. — The FUV luminosity function of our LBGs, compared with the 
LF of LBGs and QSOs from Steidel et al.(1999) and Hunt et al.(2004). Note 
that the size of error-bar of ourLF, compared to previous results, has shrunken 
significantly at bright magnitudes due to the wide area coverage of our survey. 
Dashed line is the best-fit Schechter function for z ~ 3 LBGs from Steidel 
et al.(1999), with a = -1.6, M* = -21.04 mag. Dotted line is z ~ 3 QSO 
luminosity function from Hunt et al.(2004). Solid line is the sum of the z ~ 3 
LBG luminosity function and the z ~ 3 QSO luminosity function. Our data 
points drawn with stars are consistent with the solid line. At Muv < —23.0 
mag, the QSO number density exceeds that of LBGs. The significant excess 
of LBGs at bright end over the Schechter function is most likely to be due to 
QSOs. The number density of massive (> 10"Mq) LBGs are overplotted as 
filled circles. We derived the best-fit parameters for Schechter function to be 
a = -1.6, M„ = -21.6 mag {3-dot-dashed line). 

(^1360 deg 2 ; Bentz et al. 2004). Possible candidates of these 
bright LBGs are galaxies magnified by gravitational lens (e.g., 
cB58; Williams & Lewis 1996), quasars, late-type stars, and 
low-redshift interlopers as we noted in Section 4. 1 . Among 
these possibilities, the most probable cause for bright-end ex- 
cess is QSOs at z ~ 3. According to a previous study of 
z ~ 3 QSOs (Hunt et al. 2004), the number density of faint 
QSOs exceeds that of the galaxies at the magnitude range of 
M uv < -23.0 (Figure 6). The comparison of our LBG lu- 
minosity function with the QSO luminosity function at z ~ 3 
suggests that most of the excess above Muv < -23 .0 can be 
explained by QSO. 

Only one of these Muv < -23.0 mag UV-luminous LBGs 
are detected in 24 /im image. This may be due to the depth 
of 24 fim image, and does not constrain the possibility of R- 
band bright LBGs being AGNs. If AGN-like LBGs are in- 
cluded in stellar mass analysis, the presence of AGNs may 
boost up the derived stellar masses. Still, the stellar masses 
of these bright LBGs are ~ 10 10 M Q , with the most massive 
one of 2.6 x 1O 1O M . Therefore, the excess of these "UV 
bright" LBGs does not affect the analysis of "massive", i.e., 
> 10 n M Q LBGs given in this paper. 

5.2. Properties of IRAC LBGs 

In this section, we investigate the properties of 63 IRAC 
LBGs (detected in IRAC, Sj^/xm > 6/iJy in the verification 
strip, 53.6 M m > 12/Jy in the FLS main field). We focus on 
the stellar masses and dust properties of IRAC LBGs derived 
through SED fitting. 



5.2.1. M/L ratio in rest-frame NIR 

Among 63 IRAC LBGs, 43 are identified to be more mas- 
sive than 10 n M Q . The bright rest-frame NIR flux of IRAC 
LBGs suggest that these reside at the most massive end of 
the whole LBG population, and Figure 7a shows that on av- 
erage, IRAC LBGs are at the massive end of the mass dis- 
tribution of LBGs. In Figure 7a, filled circles represent IRAC 
LBGs, filled triangles represent LBGs detected in NIR but not 
in IRAC (13), and small dots represent LBGs detected in g,R, 
and /-bands only (501). The range of the stellar masses of our 
LBGs is 2 x 10 9 M© < M* < 1O 12 M . There is no significant 
correlation between /?-band magnitudes and stellar masses. 
We estimate the uncertainty involved in stellar mass of each 
LBG to be of a factor of 2-3. 

In the local universe, the rest-frame NIR photometry can 
be a useful stellar mass indicator (e.g., Bell & de Jong 2001). 
Even for high-redshift LBGs, the same rule applies. For ex- 
ample, Rigopoulou et al.(2006) have found a correlation be- 
tween LBG stellar mass and the magnitude getting tighter as 
the wavelength of the band getting longer. In other words, 
M/L ratio of galaxies derived in the longer wavelength has 
smaller scatter compared to the value derived in the shorter 
wavelength. They used LBGs with M* < 10 n M Q , and we ex- 
tend the M/L analysis of LBGs to objects more massive than 
10 u Mq, using rest-frame NIR wavelength. Since the num- 
ber of 8 /im detected LBGs is small in our sample, we use 
4.5 /im (rest-frame 1.1 /im) flux here instead of 8.0 /xm flux 
(rest-frame TT-band). Figure 7b shows the 4.5 /im magnitude 
(the rest-frame /-band at z ~ 3) versus stellar mass. Our result 
shows that even above 10 U M Q , there is a good correlation be- 
tween the rest-frame NIR flux and the stellar mass, while no 
trend is found between the /?-band magnitude and stellar mass 
(Figure 7a). 

Despite the probable consistency of M/L ratio at rest-frame 
NIR wavelength, Shapley et al.(2005) reported that there is 
still a scatter of ~ 10 for M/L ratio of star-forming galaxies at 
Z ~ 2 when measured in rest-frame 1 .4 /zm. We also find that 
there is an order of magnitude spread in M/L ratio at 4.5 /xm. 
The variation in M/L ratio is mainly due to the differences in 
star formation history. Old galaxies have higher M/L value 
in NIR, while young galaxies show lower M/L. We exam- 
ined the M/L ratios of model galaxy templates with different 
age and star formation history, and found that M/L ratio of 
1 Gyr-old galaxy with constant star formation rate is about 
12 times smaller than 1 Gyr-old galaxy that have been evolved 
passively after a single burst (r = lOMyr). Comparing pas- 
sively evolving population of different age, the M/L ratio of 
galaxy can be more than twice higher when its age increases 
from 200 Myr to 1.5 Gyr. 

To demonstrate the above point, in Figure 7b, we mark ob- 
jects with red optical-MIR color of (R-3.6fim) > 3.5 which 
are presumably old galaxies with squares and those with high 
ongoing star formation (> 15QMq /yr) with open diamonds. 
Here, the star formation rates are derived from UV luminosity, 
corrected for the dust extinction using E(B-V) obtained with 
the SED-fitting procedure. The figure shows that actively star- 
forming galaxies have M/L ratios about 10 times smaller than 
red galaxies which have much smaller scatter in the M/L. The 
star formation activity is a major source that provides a spread 
of a factor of 10 in the mass-to-light relation even at IRAC 
wavelengths. The correlation coefficient between M4 5 Mm and 
logM* is -0.72 for all LBGs in the figure. When LBGs with 
large SFR (> \50M Q /yr) are excluded, the correlation coef- 
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ficient is -0.85. Also, when objects are restricted to galaxies 
with red (optical-MIR) colors, we find a tighter correlation 
between stellar mass and the rest-frame NIR luminosity. Sim- 
ilar effect has also been addressed in Shapley et al.(2005), 
that (R-K) > 3.5 galaxies show tighter correlation in mass- 
to-light relation among z ~ 2 star-forming galaxies. 

5.2.2. Comparison of IRAC LBGs and DRGs 

Previous studies underline the exclusive characteristics 
among the high-redshift galaxies selected by different crite- 
ria. For example, LBGs and DRGs ((J-Ks) > 2.3; Franx et 
al. 2003) are thought to represent blue and small/red and large 
systems at high redshift (van Dokkum et al. 2006). Little 
(< 10%) overlap between the two galaxy populations is ad- 
dressed (Labbe et al. 2005; van Dokkum et al. 2006). Since 
IRAC LBGs are at the massive end of the LBGs and some of 
them have red (optical-MIR) color, it is noteworthy to discuss 
whether these "IRAC LBGs" can also be selected with DRG 
selection criteria; or at least, whether their mass/color ranges 
are comparable with those of DRGs. 

As for the color, the selection cut for DRGs is (J-K) > 2.3, 
which is sensitive to the underlying old stellar population. 
There are only small number of IRAC LBGs with both J- 
and /Q-band detection, and their colors are between 0.5 < 
(J-K s ) < 1.6. With these, the IRAC LBGs cannot be selected 
by DRG criteria. However, red (optical-MIR) color cut, for 
example (R-3.6fim) > 3.5 to specify galaxies showing tight 
mass-to-light correlation, is comparable with the color cut of 
DRGs. Considering the number of (R-3.6/j,m) > 3.5 LBGs, 
we estimate that at least 20% of IRAC LBGs are candidates 
that could be selected with DRG criteria. 

Labbe et al.(2005) suggest that LBGs and DRGs are well- 
separated in (I—K s ) vs. {K s -A.5nm) diagram in Figure 1 of 
their paper. LBGs have blue (I-K s ) values (0 < (I-K s ) < 
2), and DRGs have (I-K s ) > 2, redder (K s -4.5^m) color 
than LBGs. For the IRAC LBGs with K s -band detection, we 
find that their color ranges are 1 < (I — K s ) < 4 and ((K s - 
4.5/im)) = 0.8, lying in the space between LBGs and DRGs. 



In this respect, IRAC LBGs are close to DRGs - especially in 
view of stellar masses and ages, inferred from the optical/NIR 
colors. As for the stellar mass, we show that IRAC LBGs 
have average M/L ratio of (M/Lj) ~ 0.32, while the average 
M/L ratio of DRGs is (M/L K ) ~ 0.33. Since the rest-frame 
J- and /if -band fluxes are nearly the same, the M/L ratio of 
IRAC LBGs is comparable with that of DRGs. It assures the 
similarity of IRAC LBGs and DRGs. 

The fraction of IRAC LBGs in the whole LBG population 
is ~ 6%, and it can be drawn from the known overlap be- 
tween LBGs and DRGs, ~ 10%. As is noted in Section 3.2, 
the fraction is a function of image depth - however, we see 
that IRAC LBGs represent most of "massive" high-redshift 
galaxies which could be selected from Lyman break selection 
criteria. 

5.3. Dust Properties of z ~ 3 LBGs 
5.3.1. 24 urn detected IRAC LBGs 

In this subsection, we discuss properties of twelve 24 /im 
LBGs, which are the IRAC LBGs also detected in the FLS 
MIPS image (Section 3.2). 

First, we investigate the origin of the dust emission, i.e., 
whether it comes from AGN or starburst. The combination 
of MIPS 24 /im flux and other IRAC fluxes is known to be 
a useful indicator to weed out AGNs from star-forming IR- 
bright galaxies (Egami et al. 2004; Ivison et al. 2004; Yan 
et al. 2004). Figure 8a shows the ratio of 8 /im flux to 4.5 
/im flux (Ss.o/S^.s) versus the ratio of 24 /im flux to the 8 /iin 
flux (S24/S8.0; see Figure 3 of Ivison et al. 2004). Objects 
that have large Sg.o /S4.5 and small S24 / Sa.o are considered to 
be AGNs in this diagram. Although many objects in our sam- 
ple have only upper limits in IRAC 4.5 /im, 5.8 /im or 8.0 
/zm, we examined whether the objects could be classified as 
AGN or starburst. The result shows that there are 2 possi- 
ble AGN candidates among our 24 /im LBGs, while the rest 
can be classified as star-forming galaxies that have MIR col- 
ors similar to ILLBG or cold SCUBA sources in Huang et 
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al. (2005). In a previous discussion, we mentioned that the 
bright end of the UV luminosity function could be affected 
by AGNs. Note that the two AGN-type 24 /im LBGs have 
/?-band magnitude of 23.7 mag, and 24.3 mag, therefore they 
are not directly related to the objects consisting of the bright 
end of the UV-luminosity function. 

In Figure 8b, we plot SEDs of the starburst-type and the 
AGN-type 24 /im LBGs. Solid lines are a SED of a star form- 
ing galaxy with L, R = 3 x 1O 12 L from Chary & Elbaz (2001) 
and the SED of AGN Mrk231. The top panel of Figure 8b 
confirms that most of the 24 /im LBGs have SEDs similar to 
luminous infrared galaxies with strong starburst activity. The 
rest-frame NIR part of their SEDs is nearly flat or only mod- 
erately increasing toward the longer wavelength, suggesting 
that the rest-frame NIR fluxes are mostly due to stellar light. 
The broad PAH emission features at the rest-frame 6-8 /xm 
shift to 24 /im at z ~ 3, therefore we interpret 24 /im detec- 
tion to be due to PAH emissions. On the other hand, the two 
AGN-type 24 /im LBGs have a strong power-law continuum 
typical of AGNs (Figure 8b, bottom). One of these AGN-like 
LBGs is detected in VLA 1.4 GHz with the flux of 0.59 mJy, 
which suggests that this object is a radio-loud AGN. 

In order to derive total infrared luminosities of 24 /im 
LBGs, we assumeed the SED of M82 (L IR - 3 x 10 10 L©; 
Telesco & Harper 1980) which is found to resemble SEDs 
of submm galaxies at high redshift (e.g., Lutz et al. 2005; 
MenAl'ndez-Delmestre et al. 2007) The derived infrared lu- 
minosities are of order of a few x 10 12 Lq to 10 13 Lq (Table 
3), suggesting that the 24 /im LBGs are Ultraluminous In- 
frared Galaxies. 

Alternatively, we tried to estimate L m of 24 /im LBGs 
using infrared SED templates of Chary & Elbaz (2001), in 
which the MIR flux correlates with Lj R as Lj R oc L^f l xm . With 
such a method, we obtain unreasonably high values of L IR ~ 
10 14 - 10 15 Lq for 24/mi LBGs. If 24 /im LBGs were really 
as bright as L IR ~ 10 14 - 10 15 L©, their 850 /im flux would be 
well above 10 mJy, and they shoud be detectable with submm 



observations. However, the existing submm data with a partial 
coverage of FLS show no detection of 24 /im LBGs providing 
the flux limit of < 5 mJy (Sawicki & Webb 2005; Kim et al. 
2007, in preparation). Therefore, we conclude that cautions 
are needed when deriving L IR of z ~ 3 objects using the cor- 
relation between MIR flux and IR flux such as those found in 
Chary & Elbaz (2001). 

We now investigate the dust properties of 24 /im LBGs 
by comparing SFRs derived from UV luminosity (corrected 
for dust extinction) and also from IR luminosity. When de- 
riving SFRs from UV luminosity, two methods are used to 
correct for the dust extinction. One method uses the UV 
slope P defined as f x oc X fi at 1500-2800 A (Meurer, Heck- 
man, & Calzetti 1999). Since there are no spectroscopic data 
for our LBG sample, we measured j3 from the model spec- 
tral template which is found to be the best-fit model using 
optical (u,g,R and f) photometry only. The other method 
uses the E(B-V) value derived from the SED-fitting and a 
known correlation between the UV extinction and E(B - V) 
(Aieoo = 4.39£ g (B- V); Calzetti et al. 2000). Here, our 
E(B-V) derived from stellar continuum is related with nebu- 
lar line-derived E g (B-V) as E S (B-V) = 0A4E g (B-V) (When 
we use E(B-V) value, we mean E S (B -V) hereafter). The 
conversion formula from luminosity to star formation rates 
are taken from Kennicutt (1998): 



SFR IR (M Q /yr) = 
SFR uv (M Q /yr): 



1.71xlO- lu L /R (L ) 

= 1.4 x \Q- 2i L uv {ergs/s/Hz) 



The result is presented in Table 3 and Figure 9. As is 
shown in the Table 3, the UV-derived SFRs from two meth- 
ods are consistent with each other, within a factor of 2 in 
most cases. Between the extinction-corrected UV-derived 
SFRs and IR-derived SFR, Figure 9 shows that these quan- 
tities are consistent with each other within a factor of a few, 
and with the median ratios of ~ 2 when SFRi R < 8OOM /yr 
or L IR < 5 x 1O 12 L . Note that we ignore two AGN-type ob- 
jects in this comparison. However, for galaxies with L IR > 
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TABLE 3 

Comparison of UV-derived SFR and IR-derived SFR 





object ID 


24fim flux 


UV SFR 


UV SFR C 


IR SFR 


IR luminosity 


E(B-V) 


A 1600 


comment 






(mJy) 


(M e /yr) 


(M Q /yr) 


{M^/yr) 


(£■©) 




(mag) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


c[320] FLS_ 


RJ171734.5+593548 


0.14 ±0.05 


51.7 


424.3 (1413.3) 


454.9 


2.66xl0 12 


0.36 


3.59 




c[369] FLS. 


_R_J171556.8+593833 


0.27 ±0.05 


23.8 


146.2 (163.9) 


877.2 


5.13X10 12 


0.21 


2.10 




c[421] FLS. 


_R_J171735.3+594137 


0.19 ±0.05 


35.7 


107.1 (204.8) 


617.3 


3.61 xlO 12 


0.19 


1.90 




c[583] FLS. 


.RJ171430.7+595213 


0.24 ±0.05 


16.9 


383.4 (242.4) 


779.8 


4.56X10 12 


0.29 


2.89 




c[606] FLS. 


_R_J171853.5+595325 


0.19 ±0.05 


47.6 


469.1 (473.6) 


617.3 


3.61 xlO 12 


0.25 


2.49 




c[678] FLS. 


R_J171418.8+595722 


0.56 ±0.05 


15.9 


955.4 (2723.5) 


1812.6 


1.06X10 13 


0.56 


5.59 




o[035] FLS. 


_R_J171412.0+591716 


0.61 ±0.05 


51.5 


204.5 (269.3) 


1983.6 


1.16x10" 


0.18 


1.80 


? 


o[217] FLS. 


_R_J171642.9+585733 


0.35 ±0.05 


117.8 


200.9 (155.2) 


1137.2 


6.65X10 12 


0.03 


0.30 


? 


o[275] FLS. 


_R_J171816.7+584813 


0.40 ±0.06 


16.8 


620.8 (317.6) 


1299.6 


7.60X10 12 


0.32 


3.19 




o[306] FLS. 


_R_J172045.2+585221 


1.14 ±0.06 


28.2 


35.6 (70.8) 


3710.7 


2.17X10 13 


0.10 


1.00 


AGN? a 


o[429] FLS. 


_R_J172202.1+585414 


0.45 ±0.05 


30.9 


418.7 (405.7) 


1462.1 


8.55X10 12 


0.28 


2.79 




o[504] FLS. 


_R_J171618.4+602620 


0.71 ±0.06 


15.0 


38.6 (23.8) 


2308.5 


1.35X10 13 


0.05 


0.50 


AGN? 



NOTE. — (1) ID of the object (drawn from R-band catalogs of Fadda et al.(2004)); (2) 24 fim flux in mJy; (3) uncorrected UV-derived SFR; (4) corrected UV-derived SFR using 
UV slope (3 or E{B - V) (see Section 5.3.1 for more detail); (5) IR-derived SFR; (6) IR luminosity using 24 /im flux; (7) E{B - V) derived during the SED fitting process. According to 
Calzetti et al.(2000), E S (B -V) = 0ME g (B - V), while Atem = 4.39E g (B - V) for E g (B - V) derived using nebular gas emission lines; (8) Ai 6 oo estimated from E(B - V); (9) comments 
for the objects. Objects marked as "AGN?" are suspicious to be AGNs due to their SED shapes in MIR (see Figure 8 in Section 5.3.1). Objects with "?" mark (o[035], o[217]) are 
objects with a large discrepancy between UV-derived SFR and IR-derived SFR, thus are suspected to be AGN-dominated. 

a This object is detected in VLA 1 .4 GHz image, with the flux of 0.59 mJy. 




FIG. 9. — Left: The comparison between two different extinction correction method to UV-derived SFR. Both method, using UV-slope /3 or E(B- V) are 
consistent within a factor of ~ 2. Right: The comparison between IR-derived SFR and extinction-corrected SFR. When IR-SFR is below ~ 800M Q fyr, 
extinction-corrected UV-SFR is quite consistent with IR-SFR. 



5 x 1O 12 L or SFR IR > 800M Q /yr, we find that extinction- 
corrected SFRrjyS from the both /3-correction and the E(B-V) 
methods are either systematically lower than SFR IR by a fac- 
tor of a few or have large scatter with respect to SFRi R . If 
the SFRjr values represent the true star formation rates, this 
result suggests that the extinction correction at UV tend to be 
underestimated when L IR > 5 x 10 12 L Q , consistent with the 
previous findings where it has been found that the UV-slope 
method appears to systematically underestimate the extinction 
correction in UV for the most IR luminous galaxies (Papovich 
et al. 2006; Reddy et al. 2006). 

Table 3 shows that the dust-free SFR of LBGs using 24 ^m 
flux is of order of a few hundred to a few thousand M jyr. 
Assuming that IR luminosity represents the total star forming 
activity, we derive the extinction correction at 1600 A to be 
A1600 = 2-4.5 mag for 24 fim LBGs, by comparing the SFR IR 
and SFRuv^uncorreaed- Previous works have found that the 



LBGs at z ~ 3 have a median extinction value of (Aieoo) = 1.0 
mag, distributed over A1600 = 0-4 mag (Adelberger & Stei- 
del 2000). Comparing the extinction values of 24 /im LBGs 
with the above values, we conclude that 24 [im LBGs are the 
dustiest population among LBGs. 

We also derived Ai 6 oo from E(B-V). Since we use E(B-V) 
to correct for the dust extinction of UV luminosity of LBGs 
in the following analysis, we mention here how reliable the 
extinction correction using E(B-V) would be. Mostly, there 
are ~ 2 mag scatter in Aieoo- First, the comparison between 
SFRj R and SFRuv suggest that the extinction correction using 
E(B-V) is reasonably accurate to a factor of ~ 2. Additional 
uncertainty exists in the relation between Ai o0 o and E(B-V) 
observationally (Calzetti et al. 2000) and theoretically (Witt 
& Gordon 2000). Calzetti et al.(2000) mention that the lin- 
ear relation between E(B-V) and A1600 represents the lowest 
envelope of the relation, suggesting the derived A1600 might 
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FIG. 10. — The E(B-V) versus the stellar mass. Also, overplotted on the 
figure are the data points from Rigopoulou et al.(2006). The median value 
of dust extinction in IRAC LBGs is (E(B-V)) = 0.29. On average, there 
is a correlation between dust extinction and stellar mass with a correlation 
coefficient r = 0.51. 

underestimate the dust extinction by a factor of a few. 

5.3.2. Dust property vs. stellar mass 

In Figure 10, we show a plot of dust extinction E(B-V) 
versus stellar mass. When our massive (> 1O U M ) LBGs 
are plotted together with less massive LBGs from Rigopoulou 
et al.(2006), there is a weak, but positive correlation be- 
tween logM* and E(B-V) with the correlation coefficient of 
r = 0.51. Thus, Figure 10 suggests that dustier LBGs tend 
to be more massive. It has been argued that LBGs with in- 
trinsically higher bolometric luminosity (UV+IR) are more 
massive (Reddy et al. 2006; 1.5 < z < 2.6 optically selected 
LBGs). Large amount of dust in LBG implies high infrared 
luminosity that leads to high bolometric luminosity, therefore 
our result is consistent with previous results. 

The median value for the dust extinction in IRAC LBGs 
is (E(B-V)) ~ 0.29. This is larger than (E(B-V)} ~ 0.15- 
0.20 of all LBG population (Adelberger & Steidel, 2000) and 
similar to ILLBGs of (E(B-V)} = 0.354 (Rigopoulou et al. 
2006). We already mentioned in Section 5.3.1 that the mean 
value for dust extinction is E(B-V) = 0.28 for 24/mi LBGs. 
IRAC LBGs and 24pm LBGs share the same range of dust- 
extinction, while the infrared luminosities of 24/im LBGs are 
higher than those of IRAC LBGs. By stacking the 24 fim 
images of LBGs with no individual detection in 24 rim im- 
age, we see that the average IR luminosities of IRAC LBGs 
are slightly lower, but close to those of 24/im LBGs within a 
factor of a few (presented in Section 5.3.3 in more detail). 

5.3.3. Average infrared luminosity of LBGs 

24 /^m-detected LBGs are probably the most infrared- 
luminous galaxies among all LBGs, with the total infrared 
luminosity of > a few x 1O 12 L . Since the majority of LBGs 
are not detected in 24 rim, we examined their IR properties by 
stacking MIPS 24 /im images. 

First, we examined IR properties of IRAC LBGs. We 
stacked 24 /zm images of 20 IRAC LBGs in the MIPS FLS 




FIG. 1 1 . — Left is an example of 24 (im postage image of IRAC LBG 
that is not detected in MIPS 24 image. We also present the median- 
stacked images of 20 IRAC LBGs without 24 detection in the verification 
strip (middle), and 49 IRAC LBGs without 24 detection in the main field 
(right). The size of each image is 40" X 40". The stacked image shows a 
clear signal over 8 a. 

verification strip and 49 IRAC LBGs in the main field, both 
without 24 fim detection (Figure 11). The stacked images 
show a clear existence of emission in 24 /im. When measured 
upon the stacked image, the signal to noise ratios are ~ 8 in 
both the verification and the main fields. The average fluxes 
of these LBGs, therefore, are 60.4 /tJy for 24 /tm-undetected 
IRAC LBGs in the verification strip, and 70.3 pjy for 24 fxm- 
undetected IRAC LBGs in the main field. Converting the 
expected 24 /tm flux to the total infrared luminosity assum- 
ing M82 SED, we find that the expected IR luminosity of 24 
/im undetected IRAC LBGs is about ~ 1.1 x 10 12 L© in the 
verification strip, and ~ 1.3 x 1O 12 L in the main field. The 
corresponding star formation rate is 180-220M Q /yr. The es- 
timated 24 /im flux of 60.4-70.3/zJy for IRAC LBGs is com- 
parable with the 60 /iJy cut for ILLBGs defined in Huang et 
al.(2005), who suggest the ILLBGs occupy 5% of the total 
LBGs. The fraction of our IRAC LBGs (63) in the whole 
LBGs (1088) is similar with their value. 

With the same stacking method, we also estimated the av- 
erage infrared flux of the LBGs not detected in IRAC. The 
24 ztm images of 300 LBGs in the MIPS verification field 
with no IRAC detection were stacked, showing a marginal 
detection of S/N ~ 3 or the average flux of 6 /iJy. The 
estimated infrared luminosity of these LBGs is of order of 
1.1 x 1O U L0, or the star formation rate of !9M & /yr accord- 
ing to the conversion relation. This value is consistent with 
the SFR, 10- IOOMq /yr, of typical LBGs as mentioned in 
the Introduction. At this relatively low infrared luminosity 
range, the result from Chary & Elbaz template fitting is not 
much different from the M82 scaling method; the result is 
(Li R ) = 1 .4 x 10 1 l L Q , the star formation rate of 24M Q /yr. 

5.4. Implication for Galaxy Formation 
5.4.1. Number density of massive LBGs 

The number density of massive galaxies at high redshift can 
show the straightforward evidence of early formation of mas- 
sive systems. In Figure 12, we plot the number density of 
massive (M > 1O U M ) LBGs in filled circle, and compare 
the value with other observational results (Drory et al. 2005; 
Saracco et al. 2004; Rigopoulou et al. 2006; McLure et al. 
2006) and the predictions from several hierarchical galaxy 
formation models such as semi-analytic models of Baugh et 
al.(2003; dotted line), Bower et al.(2006; dashed line), and 
the hydrodynamic simulation of Nagamine et al.(2005; large 
solid rectangle). Note that the more recent models (Bower et 
al. 2006) overcome shortcomings of the earlier model of un- 
derestimating the number density of massive galaxies at high 
redshift (e.g., Baugh et al. 2003), and are successfully pre- 
dicting the observational constraints at z ~ 2. 

We find that the number density of our LBGs with mass 
greater than 10 n M Q is $ = (1.05 ±0.15) x IQ^Mpc^. Our 
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FIG. 12. — The number density of galaxies with stellar masses > 10 u Mq 
as a function of redshift. The red filled circle shows our result for massive 
LBGs, which is quite comparable with the result of massive LBGs at the 
similar redshift range by Rigopoulou et al.(2006). Data points from other ob- 
servations are plotted with different symbols, and the prediction from recent 
semi-analytic and hydrodynamic models of galaxy formation are overplotted 
(dotted line from Baugh et al. 2003; dashed line for Bower et al. 2006; blue 
filled rectangle from Nagamine et al. 2005). 

value is consistent with the value presented in Rigopoulou et 
al. (2006), who have performed a similar study for LBGs 
over a smaller area. If we adopt ~ 20% as the fraction of 
LBGs among massive galaxies (van Dokkum et al. 2006), 
the estimated number density of all massive galaxies would 
increase to $ = (5.25 ±0.75) x 10~ 5 Mpc~ 3 , consistent with 
those from NIR-selected sample of Drory et al.(2005). QSO 
contamination at optically bright LBGs is not likely to affect 
the number density of massive LBGs since they comprise a 
small fraction of massive LBGs. 

Comparison of our result with the model predictions in Fig- 
ure 12 shows that the most up-to-date semi-analytic model 
(e.g., Bower et al. 2006) still underpredicts the number den- 
sity of massive galaxies, although the discrepancy between 
the model and the observational constraints is now much re- 
duced. Our result for the number density has uncertainty of a 
factor of a few, due to the Poisson errors and the small frac- 
tion that UV-selected galaxies comprise for the entire mas- 
sive population. However, considering that the comparison is 
made for objects at the very massive end of the mass function, 
the discrepancy could be narrowed without too much diffi- 
culty with a slight tweaking of model parameters, or by im- 
proving observational constraints with a better number statis- 
tics and measurements. 

5.4.2. Contribution of massive LBGs to the total SFR density 

In the downsizing scenario of galaxy formation, the star for- 
mation activity occurs early in massive galaxies, and late in 
less massive galaxies. Therefore, such a model predicts that 
the relative contribution of the star formation activity from 
galaxies with different masses should evolve as a function of 
redshift, with the star formation occurring more in massive 
galaxies at the higher redshift. 

In order to examine the contribution of massive LBGs to 
the total star formation rate at z ~ 3, we calculated the instan- 



taneous star formation rate using the UV luminosity function 
presented in Section 5.1, and correcting it for the dust extinc- 
tion. The FUV luminosity density is derived by integrating the 
luminosity function over the magnitude interval of the survey, 

L FUV = / L(M)$(M)dM 

JM mi „ 

while M m j„ and M„ mx indicate the minimum and maximum 
absolute FUV magnitudes. 

The estimated luminosity density is then converted to the 
star formation rate density using the conversion formula from 
Kennicutt (1998) assuming the Salpeter IMF. For the whole 
LBGs, for which we have an observational constraint on 
the faint-end slope of UV luminosity function from stud- 
ies extending to the fainter limits, we get the star forma- 
tion rate density value (before dust extinction correction) 
of p„ = 4.0 x l0~ 2 M Q yr/Mpc 3 by using M mi „ = -oo and 
M max = -10. For the star formation rate density of mas- 
sive (> 10 n M Q ) LBGs, we get a conservative estimate of 
p* = 2.3 x 10~ 4 M Q /yr/Mpc 3 by adopting M min = -23.5 and 
M max = -21.0 of our survey limit, since the faint-end slope 
of UV luminosity function of massive LBGs is not well con- 
strained. 

Correction for the dust extinction is done as the followings. 
According to the discussion in Section 5.3.1, the extinction 
parameter E(B-V) is a reasonably accurate measure of the 
dust extinction, except for IR-bright objects with very high 
star formation rate. For the SFR density from whole LBGs, 
we adopt the average dust extinction of (E(B-V)) = 0.15 
(Adelberger & Steidel 2000). This average E(B-V) value 
increases the derived SFR density by a factor of ~ 3.94 
(Calzetti et al. 2000; A(1600A) = 4.39£(B- V)), therefore the 
extinction-corrected SFR density of the whole LBGs at z ~ 3 
is p* = 1.6 X 10~ l MQ/yr/Mpc 3 . For the massive LBGs, we 
use the median E(B-V) value of massive LBGs, (E(B-V)) = 
0.29 (Figure 10), or an extinction correction factor of 14.3. 
Therefore, the corrected SFR density of massive LBGs is 
p» = 3.3 x lQ~ 3 MQ/yr/Mpc 3 . As discussed in section 5.3.1, 
the uncertainty in the extinction correction factor for 24 /im 
LBGs is about a factor of 2 to 3, but could be larger than that 
for ULIRG-type objects. We adopt the error bar of a factor of 
a few as a rough estimate of the SFRD uncertainty. 

The star formation rate density value of p, = 3.3 x 
\Q- 3 M Q /yr/Mpc 3 can be considered as a lower limit, since 
the above calculation do not include the contribution from 
massive LBGs with Muv > -21. The inclusion of the pop- 
ulation fainter than our survey limit into the calculation re- 
quires an assumption on the faint-end slope of the UV lu- 
minosity function. If we adopt a fiducial value of a = -1.6 
as the faint end slope, we get the star formation rate density 
of 1.6 x \0~ 2 M Q /yr/Mpc 3 by adopting M max = -10. This is 
likely to be an upper limit of the star formation rate of mas- 
sive galaxies, uncertain at a factor of roughly a few. Also, note 
that the derived SFR would decrease by a factor of ~ 0.25 dex 
if we adopt the Chabrier IMF (Chabrier 2003) instead of the 
Salpeter IMF (Dahlen et al. 2007). 

Figure 13 shows the SFR density of the LBGs compared 
to those of galaxies at different redshifts from other works. 
Various SFR densities in the figure are based on the UV lu- 
minosity, although sample selection methods are different. 
What is interesting in the figure is the contribution of massive 
galaxies to the total star formation rate density from z < 2 
(Juneau et al. 2005) and z ~ 3 (our result). A direct com- 
parison of our result versus Juneau et al. (2005) is possi- 
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FIG. 13. — The star formation rate density derived from LBGs in FLS. The 
open/filled circle indicates the star formation rate density from whole/massive 
(> 10 u Mq) LBGs. The points are corrected for dust extinction, using aver- 
age E(B-V) value (~ 4 and ~ 14 for whole/massive LBGs). Considering the 
contribution from LBGs below R-band magnitude limit of our study, the es- 
timated star formation rate density from massive LBGs resides in the shaded 
box. Green filled pentagons are the SFR density from galaxies whose stellar 
masses are 1O IO ' 8 M < M* < 1O 1L5 M , drawn from Juneau et al.(2005). 
Other symbols indicate the total SFR density at the corresponding redshift 
(open triangle/diamond: Dahlen et al. 2007; crosses: Schiminovich et al. 
2006; filled diamond: Hopkins 2004). Cyan filled square is the estimate of 
the SFR density at z = 3.6 using nearby massive SDSS galaxies (Panter et al. 
2006). Overplotted dotted/dashed lines are the prediction of the evolution of 
star formation rate density from all/massive galaxies along redshift (Bower et 
al. 2006). 



ble, since galaxies contributing to the most SFR in the K- 
band selected sample of Juneau et al.(2005) are star-forming 
galaxies which also appear to be UV-bright (cf. Burgarella 
et al. 2006). Despite the slight difference in stellar mass 
range (1O 1O 8 M < M* < 10 11 5 M Q in Juneau et al. 2005; 
M* > 1O U M0 in our study), our result is consistent with the 
view that the contribution of massive galaxies to the total SFR 
has steadily decreased from z ~ 3 to the present, and possibly 
providing toward an even higher contribution at z ~ 3.6 as 
suggested by Panter et al. (2006). On the other hand, a large 
uncertainty in our data point does not exclude a possible peak 
SFR activity for massive galaxies at z ~ 2 as predicted in a 
semi-analytical model of Bower et al. (2006; dashed line). 
Clearly, a better constraint is needed by extending the study 
to the fainter limit, and placing more constraints on the ex- 
tinction within high-redshift galaxies. 

Our SFR density at z ~ 3 does not include the contribution 
from the UV-faint population that cannot be selected using 
Lyman break technique, such as DRGs or submm galaxies. 
The contribution of massive galaxies to the star formation at 
z = 3 can go up even more if we include the contribution from 
the galaxies that are heavily extinguished by dust. The num- 
ber density of submm galaxies is very small, 8.9 x 10 _6 Mpc~ 3 
at 1 .8 < z < 3.6 (Tecza et al. 2004). Although there might be 
an overlap of order of 50% between LBGs and submm galax- 
ies (e.g., Chapman et al. 2005), the submm galaxies can have 
star formation rates as large as 1000 M Q /yr, therefore their 
SFR density can be almost comparable to the SFR density of 
massive LBGs. The contribution of DRGs is rather difficult 



to estimate. Forster-Schreiber et al.(2004) suggested that the 
median value of the star formation rate for DRGs in HDF- 
S is 120 M Q /yr, and their number density is about an order 
of magnitude greater than submm galaxies (van Dokkum et 
al. 2006). From the above fact, one might argue that the 
SFR density of DRG is comparable to that of massive LBGs, 
but such an argument is not valid since these two populations 
overlap with each other as mentioned in Section 5.2.2. 

In any case, our result suggests that the star formation ac- 
tivity in massive galaxies is a dominant process at z = 3 com- 
pared to the universe at z < 1, supporting the downsizing pic- 
ture of the galaxy formation. 

6. CONCLUSION 

We have selected and studied LBGs at z ~ 3 in the Spitzer 
First Look Survey area, using the multi-wavelength datasets 
consisting of deep u,g,R, and i', J, ^T s -band and the Spitzer 
IRAC and MIPS images. In total, we have found 1088 LBGs 
with R = 22.0-24.5 mag over the ~ 2.63 deg 2 area. The wide 
area coverage enables us to select a large number of bright 
LBGs with MIR fluxes which are important for studying var- 
ious properties of rare, massive LBGs. 

Particularly, we concentrated on the properties of LBGs de- 
tected in IRAC 3.6/xm over 3cr (6/iJy in verification strip, 
12/iJy in main field). These "IRAC LBGs" are on average 
massive/old/infrared sub-population of whole LBGs. Nearly 
~ 70% of IRAC LBGs are more massive than 10 n M Q . IRAC 
LBGs with the largest stellar mass have the reddest (optical- 
MIR) color, which is indicative of old stellar population. M/L 
ratio of the galaxies dominated by old stellar population is 
constant with little scatter. On the other hand, IRAC LBGs 
with large ongoing star formation increase the scatter in M/L. 

Among these IRAC LBGs, 12 LBGs were detected in 24 
/im image. The infrared luminosity of the LBGs with indi- 
vidual 24 /im detection suggests a high star formation rate of 
~ IOOOMq jyr occurring in these systems. The E(B-V) of 
LBGs, indicating the amount of dust extinction within the un- 
certainty of a factor of a few, has weak but existing correlation 
with the stellar mass. Dustier LBGs are more massive, and 
this again ensures we see the massive/infrared end of LBG 
population in the IRAC LBGs. 

With the photometric redshifts of the LBGs, the rest-frame 
UV luminosity function is constructed. The derived lumi- 
nosity function is consistent with previous studies, but with 
a much improved number statistics. We also construct the 
UV luminosity function of massive LBGs (> lO n M ), from 
which we estimate the SFR density in massive systems at high 
redshift. The star formation rate density from all LBGs at z ~ 
3, calculated in the survey magnitude interval, is />* = 1.6 x 
10~ l MQ/yr/Mpc 3 , while the star formation rate density from 
massive LBGs is p* = 3.3 x 10~ 3 - 1.6 x 10" 2 M G /yr/Mpc 3 . 
The contribution of the massive systems to the global star for- 
mation at z ~ 3 is significantly large compared to the case of 
lower redshifts. This finding suggests that the shift of star for- 
mation activity from massive systems to the smaller systems 
as the universe ages, interpreted as known "downsizing" of 
the galaxy formation and evolution. 
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